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Abstract. 



We study SU(2) Yang-Mills theory at finite temperature in the framework of the func- 
tional renormalization group. We concentrate on the effective potential for the Polyakov 
loop which serves as an order parameter for confinement. In this first status report, we 
focus on the behaviour of the effective Polyakov-loop potential at high temperatures. In ad- 
dition to the standard perturbative result, our findings provide information about the "RG 
improved" backreactions of Polyakov-loop fluctuations on the potential. We demonstrate 
that these fluctuations establish the convexity of the effective potential. 
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An understanding of strongly interacting matter at finite temperature is a promi- 
nent problem of contemporary physics that deserves to be analyzed with great ef- 
fort in view of the current and future experiments at heavy ion colliders. Since the 
forces between quarks as elementary constituents are governed by a non-Abelian 
gauge theory, already the understanding of gauge boson dynamics is an important 
challenge. In this latter case of pure gluodynamics, the expected transition to a 
deconfined phase can be studied with the aid of the Polyakov loop being the 
order parameter for this transition: 



Here t a are the generators of SU (N) and (3 denotes the inverse temperature. The 
subscript F alludes to the fundamental representation. The negative logarithm of 
the Polyakov-loop expectation value can be interpreted as the free energy of a single 
static fundamental color source 0. In this sense, an infinite free energy associated 
with confinement is indicated as (V) —>■ 0, whereas (V) ^ signals deconfinement. 

Moreover, (V) measures whether center symmetry, a discrete symmetry of Yang- 
Mills theory, is realized by the thermodynamic ensemble 0, 0] • Gauge transforma- 
tions which differ at Euclidean times Xq = and Xq = (3 by a center element of 
the gauge group change V by a phase e 2mk / N ) k integer, but leave the action and 
the functional integration measure invariant. This implies that a center-symmetric 
ground state automatically ensures (V) = 0, whereas deconfinement (V) ^ is 
related to the breaking of this symmetry. 
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1. INTRODUCTION 




(1) 



In fact, lattice simulations have not only collected strong evidence for a second 
order phase transition in SU(2) Yang-Mills theory 0, Hf, but reveal moreover 
that the critical exponents agree with those of a 3D Z(2) Ising model The 
latter corresponds exactly to the conjectured universality class obtained from the 
Polyakov-loop criterion |7[. 

In recent years, an effective theory consisting of gauge-invariant powers of the 
Polyakov loop has been developed based on mean-field arguments, see for 
an overview. Such considerations show a good agreement with lattice data 0. 
Moreover, inverse Monte-Carlo techniques have recently facilitated a precise lattice 
determination of the Polyakov-loop effective action for SU(2) [lfll ]. 

A perturbative calculation of the effective potential V for the order parameter 
was first performed by Weiss [111] . For this, it is convenient to work with the 
"Polyakov gauge" , which rotates the zeroth component Aq of the gauge field into 
the Cartan subalgebra of SU(iV); furthermore, the condition 8qA = is imposed. 
Focussing solely on the Polyakov loop, it suffices to consider a single scalar degree 
of freedom <f>(x), defined by 

A(x)^n a 6^(x). (2) 

Here n a denotes a constant unit vector in color space, e.g. n a = 5 a3 . Using this 
gauge, the order parameter is fully determined by 0, 

V& (3) 

and, consequently, the effective potential can purely be expressed in terms of 0, 
being a compact variable f3g<p e [0, 2tt] . To leading order in a derivative expansion, 
the effective potential yields [rD ] 

oo . 

/3 4 \/(/3#) = - V^cos(n/3#), (4) 

n=l 

displayed here as a Fourier-cosine series, and depicted on the left panel of Fig. 
n The potential has minima at f3gcf) = 2im and is Z(2)-symmetric, i.e., invariant 
under (3g(j) — > 2tt — j3g<p. Furthermore the order parameter is finite for /3g<f) = 2nn 
and therefore the Z(2)-symmetry is spontaneously broken and the system is in 
the deconfined phase. This perturbative result agrees with the expectation that 
perturbation theory holds for high temperature where the coupling is small. It fails 
to describe the confinement phase. 

In the confined phase, the potential should have its minimum at j3g<j) = 71 > 
implying a vanishing order parameter. A sketch of a possible form of the potential 
with a finite IR regulator (to circumvent, the convexity obstruction , see below) is 
shown on the right panel of Fig. ^for illustration. 

Various generalizations to Weiss's result have been worked out within pertur- 
bation theory, for instance, the inclusion of a magnetic background field ^2] or 
higher-order derivative expansions 0] to name a few. However, in order to inves- 
tigate the transition to the confinement phase, reliable nonperturbative tools are 



FIGURE 1. Left panel: perturbative effective Polyakov-loop potential Eq. Q [HI. A sketch 
of a possible form of the effective potential in the confining phase is given on the right panel. 



required. We will base our study on the functional (or "Exact") renormalization 
group 0] formulated in terms of a flow equation for the effective action [TBI . 

The paper is organized as follows: In Sect. Owe briefly review RG flow equations 
for Yang-Mills theories and present the flow equation for the effective potential of 
the order parameter in a propertime approximation. The flow of the potential is 
analyzed in Sect. [31 Conclusions and future directions are discussed in Sect. HI 



2. FLOW EQUATION FOR THE POLYAKOV-LOOP 

POTENTIAL 



2.1. Exact Renormalization Group 



In the flow equation approach, we consider the effective average action Tf. which 
includes all quantum fluctuations with momenta \p\ > k, with the scale k serving 
as an infrared (IR) regularization. The boundary condition for the flow equation is 
fixed at an ultraviolet (UV) scale A in terms of the bare action T\ to be quantized. 
Quantum fluctuations are successively integrated out by lowering the scale k. In the 
limit k — > 0, the full quantum effective action T^O; i-e., the generating functional 
of the 1PI Green's functions, is obtained. The flow of T^, i.e., the RG trajectory 
from the UV scale A to the deep IR, is obtained from a functional differential 
equation (IBj . 

Flow equations for gauge theories require a careful control of gauge invariance, 
because the IR regulator scale k introduces sources of gauge-symmetry breaking in 
addition to standard gauge-fixing terms. Nevertheless, standard gauge symmetry 
can be obtained in the physical limit of vanishing regulator scale k — > by control- 
ling the symmetry constraints with the aid of regulator-modified Ward-Takahashi 
identitie s |lfij |. In this paper, we employ the flow equation with the background-field 
method [17|, |l8| for a simplified control of gauge invariance within our approxima- 
tion. We work along the lines of 0, 0] where this technique has been used for 
zero-temperature gluodynamics. The flow equation for the effective average action 



reads [I?} 



Here the trace runs over all internal indices including momenta. The classical gauge 
field is given by A, whereas the background-field is denoted by A (the ghost fields 
are not displayed for brevity). T k denotes the second functional derivative with 
respect to the fluctuating fields. The regulator function R k implements the IR 
regularization at the scale k, see below. Inserting the background-field dependent 
I^ 2 ) into the regulator leads to an adjustment of the regularization to the spectral 
flow of the fluctuations as discussed in [19|, implying a potential improvement when 
it comes to approximations. 

The boundary condition for the effective action at the UV scale A consists of: 

T A [A, A] = T d [A,A} + if [A, A] + if [A, A] , (6) 

with the bare Yang- Mills action T cl . The gauge-fixing and ghost terms are given 
by 

r 9 A f [A,A] = d d x(D ll [A]a ll ) 2 , if [A 7 A] = -J d'xcD^D^Ajc, (7) 

respectively, with the quantum fluctuations defined by a = A — A. The gauge 
parameter is denoted by 

Let us briefly summarize a few properties of the regulator function R k that can 
conveniently be written as 

x 

R k {x)=xr(y), y := —— , (8) 

Z k k z 

with r(y) being a dimensionless regulator shape function of dimensionless argu- 
ment. Here Z k denotes a wave-function renormalization. Note that both R k and 
Zk are matrix- valued in field space. 

The IR regularization is implemented by the property 

lim R k (x)=Z k k 2 & r(y) v -^-. (9) 

x/k 2 — »o y 

The second and third properties of the regulator read 

lim Rdx) — and lim Rdx) — >• oo , (10) 

and ensure that the regulator is removed in the limit k — > and that the initial bare 
action is approached at the UV scale A. Moreover, the second property guarantees 
that the regulator vanishes for modes with \p\ ^> k, i.e., the theory is not affected 
by the regulator for large momenta. 



The flow equation (J3J) can be mapped onto a generalized propertime representa- 
tion, once the background field is identified with the full quantum field 1 |2lL Il9j| : 



d t r k [A=A,A) = IsTr J^f'X = \ [~ dsSTrferietp&r® ), (11) 
where the operator f(s,r)) is given by 2 

f(s,v) = Hs)(2-V) - (h(s)-g(s))(2- V ) + (H(s) -G(s))-d t . (12) 

s 

In addition, we have introduced the (matrix-valued) anomalous dimension 

r] := -d t \nZ k = -^-d t Z k . (13) 

The connection between the function f(s,rj) and the regulator function R k is given 
by the functions g(s), G(s) and H(s) which are defined as: 

%) = T^M> Hy) = J™ dsh(s)e-y s , ±H( s ) = h(s), H{0) = 0, (14) 
r(v) f°° d ~ 

^ (y) = iT^)' 9{y) = l ds ~ 9 ( s >~ ys > d: s G{s) = m ' G ^ = °- ( 15 ) 

The generalized propertime representation ([lip of the flow equation has the ad- 
vantage that the evaluation of the trace becomes considerably simplified. 



2.2. Truncated Flow Equations 

Even in the propertime form of Eq. ([ll|h the flow equation cannot be solved 
in closed form, which necessitates further approximations. For this, we truncate 
the space of action functionals down to a set of operators that are considered to 
represent the relevant degrees of freedom for the system or at least for a particular 
parametric regime of the system. 

For obtaining a first glance at finite-temperature gluodynamics, we concentrate 
on the Polyakov-loop potential, employing the simple ansatz: 

r k [A,A} = (16) 



1 This identification itself involves an approximation for in the denominator of the flow 
equation, as discussed in more detail in [ll IIJ . 

2 Terms proportional to (h(s) — g(s)) and (H(s) — G(s)) in f(s,rf) arise due to the use of the chain 
rule for d t Rk in the flow equation jSJ and manifestly represent terms arising from the spectral 
adjustment of the flow ~ dtT^ 2 \ 



Here t> M denotes the heat-bath velocity, for which we choose = S^q. In the fol- 
lowing, we neglect any running in the ghost- and gauge-fixing sectors, maintaining 
the form of Eq. (JJJ), Y 9 J' 9h = T 9 /' 9h . Furthermore, we neglect any other gauge-field 
operators except for the classical action with a wave function renormalization Z k , 
and, of course, the Polyakov-loop potential V {[v^A^) 2 ^ . In order to derive the flow 
of the potential, it suffices to evaluate the flow for a trial background field of the 
simple form 



Ad 



(n>,0) : 



(17) 



where n a denotes a constant unit vector in color space. Furthermore, we exploit 
the freedom to choose suitable wave function renormalizations in the regulator 
function for an optimal adjustment of the regulator, cf. Eq. (jHJ), 
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;is) 



for the corresponding ghost, longitudinal and transversal degrees of freedom with 
respect to the background field. In particular, we set the transversal wave-function 
renormalization equal to the background-field wave-function renormalization. The 
choice for Z^ renders the truncated flow independent of the gauge-fixing parameter 
£, so that we can implicitly choose the Landau gauge = which is known to be 
an RG fixed point 0. 

It is convenient to express the flow equation in dimensionless renormalized 
quantities, 



9k 



k^Z^f, ip = /3g<f>, v k = g 2 k~ d V k . 



(19) 



where tp 6 [0,27r]. Within our truncation, the flow equation in d— 4 dimensions 
reads 3 



perturbation theory 



+4(2- % )( I 



dxx 2 exp 



X 











1 2/ 


fk 


exp 







UJ„+X Z 



cos(nip) 



,(20) 



Here we have used the abbreviations 



Otk 



9l 
An' 



Vk 



-d t In Z k 



LOr. 



T 
zTi—n. 
k 



(21) 



Moreover, we have not displayed terms arising from oc dtT^ 2 \ e.g. terms oc dtd 2 Vk, 
on the right-hand side of the flow equation. We have dropped these terms in 



Details of the calculation will be presented in forthcoming publication pflj . 



the following preliminary numerical investigation for reasons of simplicity. As a 
consequence, the result of Eq. (J2UJ) corresponds to a standard propertime flow |2(| . 
In future work, these terms arising within the Exact RG flow will be included to 
facilitate a quantitative study of the differences between the Exact and the standard 
propertime RG in the present case. In deriving Eq. ([20)1 . we have furthermore 
employed a regulator which yields a particulary simple representation in propertime 
(or Laplace) space, 

h( s ) = 5(s-l). (22) 

At this point, it is useful to study the overlap of the present result with perturbation 
theory. In fact, we rediscover Weiss's result of Eq. (jlj) if we (i) hold the coupling 
fixed, «fc =const., (ii) set the anomalous dimension to zero, rjk = 0, and (iii) drop the 
complete second line of Eq. (J2*U|) . The resulting simplified equation is an ordinary 
differential equation, that can immediately be integrated from k = A to k = 0, 
leading us to the perturbative one-loop result of Eq. (jlj). 

Now, this observation helps estimating the nonperturbative content of the full 
Eq. (pjj): the occurrence of the running gauge coupling a fc and the fc-dependent 
anomalous dimension on the right-hand side signal the "RG improvement", i.e., 
a resummation of an infinite set of Feynman diagrams, performed by the flow 
equation. Finally, the second line of Eq. (}20|) depends on (derivatives of) the 
potential itself. This term is truly nonperturbative, being induced by fluctuations of 
the Polyakov-loop variable on top of its own potential minimum. Since fluctuations 
of the degrees of freedom associated with the order parameter become highly 
important near the phase transition, we expect that particularly terms of this 
kind can cover important aspects of the nonperturbative dynamics. 

In principle, the present truncation also allows for a rough calculation of the run- 
ning coupling gk including finite-temperature effects. For instance, the j3 function 
for the running coupling using the background field gauge reads 

d t gl = f3 gl =r) k g 2 k , (23) 

and is thus related solely to the background wave function renormalization that is 
part of the truncation. However, important features of the running coupling in the 
deep infrared require much larger truncations 1^, 24|, hence we decide to take the 



running coupling as an external input in this work. 



3. RESULTS 

Let us now discuss the flow of the Polyakov-loop potential from the perturbative 
U V regime to the IR on the basis of our minimal approximation given by Eq. (|20j) . 
In order to solve this partial differential equation, we rewrite it as an infinite set 
of coupled first-order differential equations by projecting it on a Fourier cosine 
series; this projection is naturally suggested by the fact that the Z(2) symmetry of 
the potential and its dependence on a compact variable allows for such a Fourier 




FIGURE 2. Left panel: flow of the first Fourier coefficient of V k = k d g^ 2 v k scaled by T 4 . 
Right panel: effective potential V k at different scales k for high temperature T = 5A e and g k as 
in Eq. 



expansion of the potential itself, 

v k{f) = y^ v k^ cosrup. (24) 



,(») 

n=0 



With this procedure, we obtain flow equations for the dimensionless Fourier coef- 
ficients by Fourier transforming the right-hand side (RHS) of Eq. (|20|) . 

For a numerical evaluation of the Polyakov-loop potential, we have to truncate the 
infinite set at some n = N max and set all v% with n > N max equal to zero. The 
results, shown in this work, are calculated for N max = 2 for simplicity. But we have 
confirmed that the inclusion of higher Fourier orders does not change our results 
significantly. As the boundary conditions for the flow equation at the UV scale 
A = 90 GeV, we employ the result of the one-loop calculation Eq. (jlj). 4 
Finally, our simple choice for the running coupling in this work is 



-l 

9k oc ' 1 



In (k/A c ) , (26) 



where A c denotes a "strong scale" for which we choose A c = 1 GeV for simplicity. 

In Fig. we plot the results for the flow of the coefficient Vfp (left panel) and the 
Polyakov-loop potential at different values of k as a function of p for T = 5 A c (right 
panel). The flow of the first Fourier coefficient with the running coupling of Eq. 

is compared to that for a constant coupling g k = const, and the one- loop result 
of Eq. Pjl. We find that the three different calculations are in good agreement if the 



The integration from oo to A is controlled by perturbation theory, whereas the integration from 
A to the deep IR is controlled by the nonperturbative flow. 



scale k is larger than the temperature. In this regime, the potential is built up by 
thermal fluctuations. At k ~ T, the coefficient of the nonperturbative flow-equation 
result develops a (negative) minimum and increases to zero for even smaller k. Since 
we observe a similar behaviour also for all other Fourier coefficients, the potential 
flattens for k <T, 



For the running coupling gj., this behavior is depicted on the right panel of Fig. 
El By contrast, the coefficients of the one-loop result stay finite for k — > 0. This 
flatness of the potential is, in fact, nothing but the explicit manifestation of the 
convexity property that has to hold for the effective action in general, but is missed 
by perturbation theory. 

Even though the observation of a convex potential is theoretically highly satis- 
factory, it does not tell us anything about the phase of the system. For this, the 
resulting expectation value of the Polyakov loop in the IR is relevant. At this point, 
we should stress that, even for a finally flat potential, the expectation value does 
not remain undetermined, but is well-defined by its k — > + limit (the effective 
potential is not flat for any nonzero k) . However, as we can read off from Fig. 
we do not observe a sign change of the Fourier coefficients under the flow; hence 
the minima of the effective potential are Z(2) symmetry breaking and thus corre- 
spond to the deconfinement phase. In the present simple truncation, this holds true 
even for lower temperatures. Consequently, our truncation so far is only capable of 
describing deconfined dynamics. 



We have presented first steps towards a nonperturbative study of the Polyakov- 
loop potential based on an RG flow equation. Our intention so far mainly was to 
demonstrate the capability of our approach as a matter of principle by choosing 
a minimalistic approximation scheme. Already at this level, we observe important 
nonperturbative features, such as the backreactions of order-parameter fluctua- 
tions on top of the potential minimum and convexity of the effective potential. 
Confronting our first results with phenomena, we find that our truncation oversim- 
plified the system, since we find deconfined dynamics on all scales. With hindsight, 
this is not too surprising, since, in addition to the Polyakov-loop dynamics, the 
gluon (and ghost) sector is basically approximated by its classical form, retaining 
"perturbative" gluon degrees of freedom all the way down to k = 0. From this point 
of view, our present minimalistic truncation leads to self-consistent results. 

Future work will extend the present approach in various directions: on a technical 
level, the influence of the (so far neglected) terms from spectral adjustment oc 

(2) 

dtT v k has to be studied. This is an interesting task by its own, since it allows 
a quantitative analysis of the importance of these terms, distinguishing standard 
propertime flows from (background-field approximated) Exact RG flows. On a 
more conceptual level, the truncation has to be extended to include a larger 




and 




(27) 



4. CONCLUSIONS AND OUTLOOK 



class of gluonic (and/or ghost) degrees of freedom, covering important aspects 
of nonperturbative dynamics and allowing for an interplay with the Polyakov-loop 
sector. For instance, a gluonic potential Wki^F^F^) replacing the present simple 
ansatz oc F^F^ leaves room for the description of a nontrivial magnetic sector 
including gluon condensation as well as complex dynamics in combination with 
the Polyakov loop. Future work in this direction is in progress. 
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